Volksbank

Load Libraries

library(geojsonio)
library(leaflet)
library(magrittr)
library(htmlwidgets)
library(RColorBrewer)
library(tidyverse)
library(rgdal)

Amsterdam

Load data

# Read shapefile: Spatial Polygon DB
amsterdamData <- geojson_read("geojson.json", method= "local", what = "sp")

Create maps

# Show Amsterdam building projects on a Leaflet map
# fasepal <- colorFactor(c("DeepSkyBlue", "MediumSpringGreen ", "Orange", "Tomato"), amsterdamData$Fase)

fasepal <- colorNumeric(c("Tomato", "Orange", "MediumSpringGreen", "Skyblue"), amsterdamData$Start_bouw)

leaflet(amsterdamData) %>%
  setView(lng = 4.8945, lat = 52.3667, zoom = 12) %>%
  addProviderTiles("nlmaps.standaard") %>% 
  addProviderTiles("Esri.WorldGrayCanvas") %>%
  addPolygons(stroke = TRUE, color = "white", weight="1", smoothFactor = 0.3, fillOpacity = 0.7, fillColor = ~fasepal(Start_bouw), popup = ~Projectnaam) %>%
  addLegend(pal = fasepal, values = ~Start_bouw, opacity = 0.8, title = "Start bouw", position = "topright", labFormat = labelFormat(big.mark = ""))

Load CBS data

# Alleen eerste keer
# fileName <- "buurten.geojson"
# wijkenenbuurten <- geojson_read(fileName, what = "sp")
# wijkenenbuurten <- wijkenenbuurten[na.omit(wijkenenbuurten@data$GM_NAAM == "Amsterdam"),]
# writeOGR(wijkenenbuurten, "amsterdamBuurten.geojson", "GeoJSON", driver="GeoJSON")

amsterdamBuurtenData <- geojson_read("amsterdamBuurten.geojson", what = "sp")
amsterdamBuurtenData65 <- amsterdamBuurtenData[amsterdamBuurtenData$P_65_EO_JR >0 ,]

agepal <- colorNumeric(c("Tomato", "Orange", "MediumSpringGreen", "Skyblue"), amsterdamBuurtenData65@data$P_65_EO_JR)

leaflet(amsterdamBuurtenData65) %>%
  setView(lng = 4.8945, lat = 52.3667, zoom = 12) %>%
  addProviderTiles("nlmaps.standaard") %>% 
  addProviderTiles("Esri.WorldGrayCanvas") %>%
  addPolygons(stroke = TRUE, color = "white", weight="1", smoothFactor = 0.3, fillOpacity = 0.7, fillColor = ~agepal(P_65_EO_JR), popup = ~BU_NAAM) %>%
  addLegend(pal = agepal, values = ~P_65_EO_JR, opacity = 0.8, title = "% 65 jaar en ouder", position = "topright") 

amsterdamBuurtenData <- geojson_read("amsterdamBuurten.geojson", what = "sp")
amsterdamBuurtenDataWOZ <- amsterdamBuurtenData[amsterdamBuurtenData$WOZ > 0 ,]

wozpal <- colorNumeric(c("Tomato", "Orange", "MediumSpringGreen", "Skyblue"), amsterdamBuurtenDataWOZ$WOZ)

leaflet(amsterdamBuurtenDataWOZ) %>%
  setView(lng = 4.8945, lat = 52.3667, zoom = 12) %>%
  addProviderTiles("nlmaps.standaard") %>% 
  addProviderTiles("Esri.WorldGrayCanvas") %>%
  addPolygons(stroke = TRUE, color = "white", weight="1", smoothFactor = 0.3, fillOpacity = 0.7, fillColor = ~wozpal(WOZ), popup = ~BU_NAAM) %>%
  addLegend(pal = wozpal, values = ~WOZ, opacity = 0.8, title = "WOZ-waarde (x 1.000 EUR)", position = "topright") 

Combine maps

leaflet() %>%
  setView(lng = 4.8945, lat = 52.3667, zoom = 12) %>%
  addProviderTiles("nlmaps.standaard") %>% 
  addProviderTiles("Esri.WorldGrayCanvas") %>%
  addPolygons(data = amsterdamBuurtenDataWOZ, stroke = TRUE, color = "white", weight="1", smoothFactor = 0.3, fillOpacity = 0.4, fillColor = ~wozpal(WOZ), popup = ~BU_NAAM) %>%
  addLegend(data = amsterdamBuurtenDataWOZ, pal = wozpal, values = ~WOZ, opacity = 0.8, title = "WOZ-waarde (x 1.000 EUR)", position = "topright") %>%
  addPolygons(data = amsterdamData, stroke = TRUE, color = "white", weight="1", smoothFactor = 0.3, fillOpacity = 0.9, fillColor = "DodgerBlue ")

leaflet() %>%
  setView(lng = 4.8945, lat = 52.3667, zoom = 12) %>%
  addProviderTiles("nlmaps.standaard") %>% 
  addProviderTiles("Esri.WorldGrayCanvas") %>%
  addPolygons(data = amsterdamBuurtenData65, stroke = TRUE, color = "white", weight="1", smoothFactor = 0.3, fillOpacity = 0.4, fillColor = ~agepal(P_65_EO_JR), popup = ~BU_NAAM) %>%
  addLegend(data = amsterdamBuurtenData65, pal = agepal, values = ~P_65_EO_JR, opacity = 0.8, title = "% 65 jaar en ouder", position = "topright") %>%
  addPolygons(data = amsterdamData, stroke = TRUE, color = "white", weight="1", smoothFactor = 0.3, fillOpacity = 0.7, fillColor = "DodgerBlue ")

Oss

# Alleen eerste keer
# fileName <- "buurten.geojson"
# wijkenenbuurtenOss <- geojson_read(fileName, what = "sp")
# wijkenenbuurtenOss <- wijkenenbuurtenOss[na.omit(wijkenenbuurtenOss@data$GM_NAAM == "Oss"),]
# writeOGR(wijkenenbuurtenOss, "ossBuurten.geojson", "GeoJSON", driver="GeoJSON")

wijkenenbuurtenOss <- geojson_read("ossBuurten.geojson", what = "sp")
wijkenenbuurtenOss <- wijkenenbuurtenOss[wijkenenbuurtenOss$WOZ > 0 ,]

leaflet(wijkenenbuurtenOss) %>%
  setView(lng = 5.5140, lat = 51.7612, zoom = 12) %>%
  addProviderTiles("nlmaps.standaard") %>% 
  addProviderTiles("Esri.WorldGrayCanvas") %>%
  addPolygons(stroke = TRUE, color = "white", weight="1", smoothFactor = 0.3, fillOpacity = 0.7, fillColor = ~wozpal(WOZ), popup = ~BU_NAAM) %>%
  addLegend(pal = wozpal, values = ~WOZ, opacity = 0.8, title = "% 65 jaar en ouder", position = "topright")